Evidence for Dirac Fermions in a honeycomb lattice based on silicon 
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Silicene, a sheet of silicon atoms in a honeycomb lattice, was proposed to be a new Dirac-type 
electron system similar as graphene. We performed scanning tunneling microscopy and spectroscopy 
studies on the atomic and electronic properties of silicene on Ag(lll). An unexpected y/Z x \/3 
reconstruction was found, which is explained by an extra-buckling model. Pronounced quasi-particle 
interferences (QPI) patterns , originating from both the intervalley and intravalley scattering, were 
observed. From the QPI patterns we derived a linear energy-momentum dispersion and a large 
Fermi velocity, which prove the existence of Dirac Fermions in silicene. 

PACS numbers: 68.37.Ef, 73.22.-f, 61.48.-c, 71.20.Mq 



Group IV (Si, Ge) analogs of graphite have been dis- 
cussed for a long history even before the synthesis of 
isolated graphene 0, and recently there have been re- 
nev^fed interest in this topic due to the novel concepts 
and applications brought by graphene. The silicon ver- 
sion of graphene in which Si atoms replace C atoms in 
a two-dimensional honeycomb lattice is named silicene 
[H, ■ Theoretical calculations show that silicene has 
also graphene-like electronic band structure, support- 
ing charge carriers behaving as massless Dirac Fermions 
[1, [^. Compared with graphene [1, 0], silicene has a 
larger spin-orbit coupling strength, which may lead to 
larger energy gap at the Dirac point and favor detectable 
quantum spin Hall effect (QSHE)[5l, [13, Currently 



QSHE has been realized only in HgTe-CdTe quantum 
wells and further studies have been hindered by the chal- 
lenging material preparation [l^ . The easy prepara- 
tion and compatibility with silicon-based nanotechnol- 
ogy makes silicene particularly interesting for applica- 
tions like QSHE devices. 

Despite the rapidly increasing amount of theoretical 
works on silicene, there have been only a few experiments 
on silicene or silicene nanoribbon [H, 0, 0] ■ Monolayer sil- 
icene film has been successfully grown on Ag(lll) [4] and 
ZrB2 substrates [l3| . and scanning tunneling microscopy 
(STM) study revealed hexagonal honeycomb structure, 
which was distinct from known surface structures of bulk 
silicon, and resembles that of graphene Q. However, 
the reported Si-Si distance of 17% shorter than that for 
the bulk and for the theoretical model implies unrealistic 
high compression of the silicene lattice, which remains to 
be confirmed and understood [4]. Apart from the prepa- 
ration and structural studies, there is still no experiments 



on the electronic structure of monolayer silicene. Experi- 
mentally establishing a common understanding of the ba- 
sic atomic and electronic properties of silicene is therefore 
highly desirable. 

In this Letter we report a study on silicene by low 
temperature STM and scanning tunneling spectroscopy 
(STS). We observed an unexpected y/2> x -s/S reconstruc- 
tion on silicene surface, in contrast to the 1x1 structure 
reported previously Despite the structural difference, 
the electronic property measured by STS is consistent 
with theory very well. For examples, quasi-particle in- 
terference (QPI) patterns suggesting intervalley and in- 
travalley scattering of charge carriers were observed, and 
a linear energy-momentum dispersion relation and a large 
Fermi velocity were derived. Our results provide a solid 
basis for further studies on electronic property and device 
applications of silicene. 

Experiments were performed in a home-built low- 
temperature STM equipped with in-situ sample prepara- 
tion facilities. A clean Ag(lll) surface was prepared by 
cycles of argon ion sputtering and annealing. Silicon was 
evaporated from a heated Si wafer with a deposition flux 
of about 0.05 ML/min, and the Ag(lll) substrate main- 
tained at wSOO K. The differential conductance (dl/dV) 
spectra were measured as the in-plane ac component in 
the tunneling current with a lock-in amplifier by superim- 
posing an ac voltage of 10 mV and 676 Hz on the given 
dc bias of the substrate-tip gap. We performed STM 
measurements at both room temperature and 77 K, and 
found no difference. All the data presented in this paper 
were taken at 77 K. 

The growth of Si on Ag(lll) is very sensitive to the 
Si coverage and the substrate temperature . We have 
explored a wide range of Si coverage and substrate tem- 
perature, and several metastable phases were found at 
low temperature range 15|. On the other hand, at sufR- 
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ciently high substrate temperature and small Si coverage, 
silicon routinely forms one-atom thick, closely packed is- 
lands which we identify as silicene sheets. Typical STM 
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FIG. 1: (color online) The STM image of a large area (65 nm 
65 nm) consisting of a sheet of silicene on Ag(lll) crossing two 
substrate steps, (b) The line profile as indicated by the black 
line in (a) shows that the island is of one atom thick, (c) The 
high-resolution STM image (10 nmxlO nm) of the silicene 
surface taken at tip bias 1.0 V. The honeycomb structure is 
clearly observed, (d) The line profile as indicated by the black 
line in (c) showing both the lateral and vertical corrugation 
of the structure observed by STM. 
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FIG. 2: (color online) (a), (b) The top view and side 
view of the lattice geometry of the y/S x y/S superstructure 
(ABA) and low-buckled phase (AB), respectively. Note that 
in the \/3 x \/3 superstructure a Si atom with planar coor- 
dinate (2/3, 2/3) is pulled downward. The red, yellow and 
green balls represent the A, B and A Si atoms, respectively. 
(c)The structural phase transition diagram of silicene depend- 
ing on the lattice constants of y/S x y/S superstructure, (d) A 
larger schematic model illuminating the honeycomb structure 
of y/S X y/S reconstructed silicene. 



topographic image of the thus obtained monolayer islands 
is shown in Fig. 1(a). The line profile across the edge of 
a silicene island indicates that it has a height of 2.63 A, 
corresponding to one atom thickness. Such islands can 
run across steps of the Ag(lll) substrate without losing 
the continuity of the atomic lattice, similar to that of a 
graphene sheet falling on a stepped surface. As expected, 
the high resolution STM image in Fig. 1(c) shows a hon- 
eycomb structure reflecting the three-fold symmetry of 
this film. However, surprisingly, after a very careful cal- 
ibration we derived a periodic constant of 0.64±0.01 nm 
(shown in line profile in Fig. 1(d)), which is approximately 
y/3 a (a is the lattice constant of silicene, ~0.38 nm, pro- 
posed in several theoretical works d, Q). Such a large 
difference is impossible a result of strain-induced lattice 
expansion of compression. It is most likely that we ob- 
served a (\/3x a/3)R.30° superstructure, which frequently 
shows up in three-fold symmetric systems. 

To account for this observation, we have carried 
out first-principles calculations using the projector 
augmented wave (PAW) pseudopotential method and 
Perdew-Burke-Ernzerhof (PBE) exchange-correlation 



potential [16| implemented in the VASP package jl7|. 
The convergence criteria for energy and force were set 
to 10~^eV and O.OOleV/A, respectively. We performed 
an extensive search for different geometric configurations 
with {y/3 X -y/3)R.30° periodicity, and concluded that a 
{y/3 X \/3)R30° superstructure cannot be stabilized in a 



free-standing, fully relaxed silicene model. We also con- 
sidered the incorporation of Ag atoms in -\/3 x -\/3 sites 
in the silicene lattice, but it was proven to be energeti- 
cally unfavorable. The period constant 0.64 nm cannot 
also be resulting from the commensuration between the 
lattices of silicene and Ag(lll). Actually the theoreti- 
cal modeling only gives a y/7 x -s/? superstructure (with 
respect to the 1x1 Ag substrate). We also exclude the 
possibility that the honeycomb (y/Sx a/3)R30° structure 
is simply an intervalley QPI pattern which was known 
in graphene to create a (y/S x •\/3)R.30° pattern [3]. In 
fact, we observed such interference pattern (to be shown 
in the following context), but they can be distinguished 
from the (y/S x -\/3)R.30° superstructure observed here, 
indicating that they have different origins. 

Finally, we found a clue by comparing the obtained 
periodicity, 0.64 nm, with the theoretically proposed val- 
ues, all around 0.67 nm 0, [l^. The experimental 
value is about 4% smaller than the theoretical ones, in- 
dicating that the substrate may exert some influence to 
result in the contraction of the silicene lattice. Although 
the detail mechanism for such contraction still remains 
to be understood, for a phenomenological model we can 
fix the periodicity constant in a series of contracted sil- 
icene, and fully relax the inner coordinates of Si atoms. 
Interestingly, we found that the (VS x V3)R30° recon- 
struction can indeed be stabilized in a contracted lattice. 
The Si atoms in previous theoretical low-buckled silicene 
(Fig. 2(a)) have two different heights, which we named 



AB configurations. To understand the experimentally 
observed (-s/S x VS) superstructure, we considered a con- 
figuration with extra buckling of the Si atom with planar 
coordinate (2/3, 2/3) downward in one (v^x V^) silicene 
unit cell, which is named ABA configuration, as shown 
in Fig. 2(b). We carried on a comprehensive study of the 
structure and stability of a series of contracted silicene 
structures from the AB and ABA configuration with and 
without small perturbation as the initial structures, re- 
spectively. Fig. 2(c) shows the structure phase transition 
diagram of the (-^3 x ^/3) superstructure and the low- 
buckled AB phase with minimum energy and stability, 
through structural optimization by keeping the periodic- 
ity constant in the contracted region. When the lattice 
constant is smaller than 6.25 A, both the ABA configu- 
ration and the AB configuration as the initial structures 
return to the (-y/S x -\/3) superstructure (ABA configu- 
ration). In contrast, when the lattice constant is larger 
than 6.45 A, an initial ABA configuration will sponta- 
neously transfer to the low-buckled AB configuration. 
The \/3 X \/3 superstructure (ABA configuration) is ro- 
bust and stable when the lattice constant is smaller than 
6.35 A. Taking into account the scale error, our exper- 
imental silicene sample may be in the ABA phase with 
the (-v/S X ^/3) superstructure. 

Silicene has been predicted to have Dirac-type electron 
structure around the Fermi energy [^, Q , similar to that 
of graphene. This however remains to be experimentally 
confirmed. On the other hand, measuring the charac- 
teristic electronic structure of our silicene film can also 
provide an additional proof of the basic atomic structure 
of our film. For graphene, the Dirac cones are located at 
high symmetric K points in the 2D Brillouin zone. The 
constant energy contours in reciprocal space cut through 
the electron or hole conical sheets near Ep resulting in 
small circles centered at the K points (shown in Fig. 
3(d)). Free carriers can be scattered within the small 
circles (intravalley scattering) or between circles (inter- 



valley scattering) 23-22|, resulting in QPI patterns in 
real space, which have been observed STS. Similar QPI 
patterns are expected to exist for silicene based on the 
same analysis. 

We have observed short-wavelength interference pat- 
terns in STM images of silicene, consistent with an in- 
tervalley scattering. From Fig. 3(d) one can see that 
the wave vectors of intervalley scattering, q2 , is close 
in length to the wave factor K]^(2)j which has a length 
of l/\/3 in reciprocal space, corresponding to a (V3 x 
\/3)R30° periodicity in real space. Sucha(y3xV3)R30° 
interference pattern had been observed in areas near step 
edges and defects which serve as scattering centers, as il- 
lustrated in Fig. 3(a)-(c). The red circle in (b) is drawn 
near the step edge, and in (c) it is along the close-packing 
direction away from the point defect. Within the red cir- 
cles the surface is imaged as close-packed protrusions. 
In contrast, within the white circles the surface is im- 




FIG. 3: (color online) (a) The STM image (50 nm x 45 nm) of 
a silicene island taken at tip bias -1.1 V. The white boxes label 
two typical defect sites, namely step edge and point defect, 
where scattering patterns with hexagonal close packing struc- 
tures were observed, as enlarged in (b) and (c). (d) Schematic 
of 2D Brillouin Zone (Blue lines), constant energy contours 
(grey rings) at K points, and the two types of scattering vec- 
tors: qi (intravalley, short red arrow) and q2 (intervalley, long 
red arrow). 



aged as honeycomb structure, as commonly observed in 
large area silicene islands. Therefore, we can find a clear 
phase shift between the atomic corrugation and the QPI 
pattern, although both of them exhibit x •\/3)R.30° 
periodicity. The QPI pattern originating from the scat- 
tering center extends for a length of only a few nm. Such 
a characteristic decaying length cannot make us assign 
the (-s/S X •\/3)R30° structure observed on larger islands 
simply as an QPI pattern. The observation of QPI pat- 
terns, consistent with analysis based on the theoretical 
band structure of silicene is an additional proof that the 
underlying atomic structure of our film is graphene-like, 
with only some buckling that does not change the basic 
electronic structure of the film. 

To investigate the electronic structure of silicene in 
more detail, we performed STS measurements (dl/dV 
curves and maps) on the film. Typical dl/dV curve taken 
at 77 K is shown in Fig. 4(a). Beside the pronounced peak 
at 0.9 V, there is a small dip located at about 0.5 V, 
which is attributed to the position of Dirac point (DP) 
of silicene. The dip in dl/dV curve corresponding to 
DP is not much obvious compared with that of graphene 
24l |. which is due to the pronounced electronic den- 



23 



sity of states (DOS) of the underlying Ag(lll) substrate 
superimposed in the dl/dV spectra. The Fig. 4(c)-(e) are 
dl/dV maps of a silicene island (second layer of silicene) 
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(1.2±0.1) X 10^ m/s. The k=0 energy intercept gives the 
Dirac energy, 'Ep-'Eu = 0.52±0.02 eV, in consistent with 
the position of DP in dl/dV spectra (Fig. 4(a)) very 
weU. The hnear E-k dispersion proves the existence of 
Dirac cone in electronic band structures of sihcene. The 
surprising lar ge F ermi velocity, comparable with that of 
graphene 2^, 2^, suggests the prospective applications 
comparable with those have been proposed or realized in 
graphene. 

In summary, we have obtained silicene film on Ag(lll) 
and observed an unexpected a/Sx -^3 buckled superstruc- 
ture. QPI patterns resulting from both intervalley and 
intravalley scattering were found, and the study of the 
linear energy-momentum dispersion and the measured 
Fermi velocity, as high as 10^ m/s, proves that quasi- 
particles in silicene behave as massless Dirac fermions. 
Based on these results, many further interesting stud- 
ies may be proceed, including gap opening in silicene, 
which can realize the quantum spin Hall effect [HI, and 
high temperature electron-phonon superconductivity in 
hydrogenated silicene [26. 27|. 
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FIG. 4: (color online) (a) dl/dV curves taken at 77 K. The 
position of Dirac point (DP) is labeled, (b) The STM image 
(40 nm X 40 nm) of IML silicene surface containing an island 
of second layer taken at tip bias -1.0 V. (c), (d) and (e) dl/dV 
maps of the same area as (b) taken at tip bias -0.8 V, -0.7 V 
and -0.6 V, respectively, (f) Energy dispersion as a function 
of K for silicene determined from wavelength of QPI patterns. 
The inset shows an schematic drawing of the overall band 
structure, with the relative location of DP, Ef and our data 
points (red line) 



on the single layer of silicene, reflecting the distribution 
of LDOS in real space. Wave-like QPI patterns near the 
boundary of the island were clearly observed. The wave- 
length changes as a function of tip bias voltage. As the 
bias increases from -0.4 V to -1.1 V, the wavelength de- 
creases correspondingly from 2.8 nm to 1.6 nm. 

The 2D constant-energy contours in reciprocal space 
(Fig. 3(d)) are used to understand the QPI patterns. The 
wave factor qi of intravalley scattering connect points 
within a single constant-energy circle and determine the 
observed long wavelength interference pattern. In order 
to deduce the quasipartical energy-momentum dispersion 
relation, we drew E(k) curve in Fig. 4(f), where k is 
the radius of constant-energy circle at K point with 2k— 
qi. The values of qi are determined by measuring the 
wave length of QPI patterns in dl/dV maps. We found 
AC varied linearly with energy, with Fermi velocity Vp = 
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